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Abstract: We study the structure of the phase space in Hofava-Lifshitz theory. With the 
constraints derived from the action, the phase space is described by five fields, thus there is 
a lack of canonical structure. The Poisson brackets of the Hamiltonian density do not form 
a closed structure, resulting in many new constraints. Taking these new constraints into 
account, it appears that there is no degree of freedom left, or the phase space is reduced 
to one with an odd number of fields. 
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1. Introduction 

Hofava recently proposed a gravity theory with asymmetry between time and space [l], [2j. 
This theory is non-relativistic in the UV limit, thus it is hoped that it is UV finite. It is 
similar to a scalar field theory of Lifshitz || in which the time dimension has weight 3 if 
a space dimension has weight 1, thus this theory is called Hofava-Lifshitz gravity. Hofava 
argued that it is superficially renormalizable based on power counting and may flow to 
Einstein's general relativity in the IR region. This work has stirred up a surge of research 
on possible applications of the theory to cosmology and black hole physics. We have no 
intention to be complete in offering the literature, for those interested in a list of these 
papers, we refer to a most recent paper ||. 

When any theory claims to be a renormalizable field theory of gravity, one must exercise 
great care, especially with a theory without general covariance to begin with, since we know 
that general covariance is hard to avoid in a theory with a massless spin 2 particle. One may 
already poses questions at face value. If a theory can flow to Einstein's theory only when 
a cosmological constant is introduced, how can one avoid this cosmological constant? 1 If 
one fine-tunes parameter A in the kinetic term in the Hofava action to 1/3 thus makes the 
cosmological constant vanish, then how A can flow to 1 in Einstein theory? And, how a field 
theory of gravity can explain the fact that the maximal entropy of a region is proportional 
to the area of the surface surrounding it? 

With the above questions in mind, we begin a study on constraints in this theory. 
One easily comes to doubt that whether the system of constraints of the theory makes 
sense, since not all constraints correspond to local symmetries. Hofava retains all diffeo- 
morphism symmetries in space, but gives up on time local symmetry. In a constrained 
system, constraints are normally generators of symmetries, thus they are guaranteed to 
form a closed system under the Poisson bracket. Now, the constraints corresponding to 
the lapse function have no corresponding symmetries, it is natural that they will generate 
new constraints under the Poisson bracket. Indeed they do, as we will show shortly. 

One may simply choose not to impose these new constraints. If so, then there are four 
families of constraints: one corresponding to the lapse function N, three corresponding 

1 Some papers also pointed out the problem related to cosmological constant M, m. 
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to shift functions N l . In the ADM canonical formalism (especially suitable for the non- 
relativistic theory of Hofava), there are 12 fields on the phase space, 6 are the spatial 
metric components gij, 6 are their canonical momenta. Upon imposing four constraints 
and three gauge symmetries, the phase space is described by five functions, and there is no 
symplectic structure on this phase space. 

Thus we need to impose new constraints derived from the Poisson brackets, but we 
shall see that there are too many new constraints, thus it appears that no degree of freedom 
is left, or the phase space is reduced to a smaller one still described by an odd number of 
fields, this seems to us a fatal problem in Hofava's theory. 



2. Non-closure of constraints algebra and New constraints 

We start with the relativistic metric g^ v in the usual ADM decomposition, 

-N 2 + NiN 1 Ni 



(2.1) 



for any theory of metric to be a theory of gravity, the lapse function N must be a function of 
both space and time, since Newtonian potential is embedded in it. Even though Newtonian 
potential can be included in the shift functions by choosing certain special gauge, to obtain 
Newtonian equation which determines Newtonian potential, the lapse function iV must be 
a function of both space and time. Without Newtonian equation as a local constraint, 
many unphysical solutions will emerge. In UV region the action in Hofava-Lifshitz theory 
takes the following form 

S = ^ J dtcPx^gNUKijK** - \K 2 ) - ^-C^}, (2.2) 

where for simplicity we consider the first action proposed in ||] in this section, we postpone 
a discussion of the second action (which may be generated by flowing the above action 
and flow to Einstein action in the IR region) in |2| to sect. 3, where similar calculations are 
carried out. The first two terms in the above action comprise the kinetic term and the 
last term is the potential term. This action satisfies the detailed balance principle, since 
C ij = e ikl V k (Rj - \R8{) is the Cotton tensor and can be obtained from the variation of a 
3-dimensional action || 

5W 

y/gV = (2.3) 

where W[<7ij] is the 3-dimensional Chern-Simons action. Kij = ^{(jij — VjiVj — VjiVj) are 
components of the extrinsic curvature of the 3-dimensional hyper-surface with constant t 
and K = g 13 Kij is its trace. k\y is a constant in this theory with mass dimension 1. 

In the following, we focus on the case of A = 1 so that the kinetic term in the action 
looks like that appearing in Einstein-Hilbert action. Our calculation demonstrates that 
different value of A will not change the final results qualitatively, so we can make this 
choice. Hereafter, we will adopt the units in which 16irG = 1, for convenience. In the 
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Hamiltonian formalism, we first compute the conjugate momenta of N, Ni and gij denoted 
by 7r,7r l and it 13 , respectively. They have the explicit form 

vr = *l = 0, (2.4) 
^ = 6S_ = 

^ = ^- = V9(K ij -g ij K). (2.6) 

The basic Poisson brackets of the canonical variables are 

{N(x), vr(y)} Pb = 5 3 (x-y), (2.7) 
{Nt(x), ^'(y)} Pb = <^ 3 (x-y), (2.8) 



{ 9ij (x),7r kl (y)} Ph = (Sf5 l j + 5lS^(x-y), (2.9) 



where we have adopted the convention of Poisson brackets between canonical variables 
used by (4j. One should be aware of the fact that ix^ is not a tensor under coordinate 
transformation, but behaves as a tensor density. In other words, Tiijj ' J~g is a tensor. Then 
after performing a Legendre transformation, the Hamiltonian can be derived as 

H = j TT ij gijd 3 x - L = J d 3 x(NH + N l Hi), (2.10) 

with 7i and Hi given by 

n = ^p^Gijki^ + -j^-VgCijC**, (2.ii) 

Hi = -Igud^i - (2d kgij - d i9jk )^ k . (2.12) 

where Q ijkl = \{g ik gji + gugjk - gijgu) is the inverse of de Witt metric. Since Eqs.(j2l|) 
and ( |2.5| ) show that the canonical momenta of N and ./V* always vanish, the two equations 
below are always satisfied 

ir = {n, H} Ph = H, (2.13) 

V = K, H} Ph = Hi, (2.14) 

provided the following constraints 

H = 0, (2.15) 

Hi = 0. (2.16) 

Above equations are usually called as the Hamiltonian constraint and momentum con- 
straint, respectively. Once the Poisson bracket is added into the system as a part of 
structure, consistency requires that the Poisson brackets among the constraints only gen- 
erate constraints. As discussed in the introduction, we anticipate that new constraints 
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are generated by the Poisson brackets among Hamiltonian constraints. To compute the 
Poisson brackets of each pair of constraints, it is helpful to first compute 

{ J d 3 x'c k n k , gij } Ph = -c k d k9iJ - 9jk diC k - gik d^ h , (2.17) 

{ J d^'( k H k ,^} Ph = -d k (^C k )+7T^d k C + 7T lk d k C^ (2.18) 

which reveal the role of Hi's as generators of 3-dimensional coordinate transformation, and 
{ J d 3 x.'r]H, gij } Ph = -^=(2g u g jk - gij g k i)ir kl = -2rji%. (2.19) 
Then one obtains the following two Poisson brackets straightforwardly, 

{ j tfxGHi, J <rycpt,}v-, = J d 3 x(CA(l - (idi(i)n k , (2.20) 

{ J d^CHi, J d 3 y V H} Ph = J d 3 xCd tV n, (2.21) 



Eq.( 2.2(3j ) tells us that Hi's form a Lie algebra corresponding to 3-dimensional diffeomor- 



phism group. Eq. Q2.21 ) is a reflection of the fact that H behaves as a scalar density under 



this transformation. Actually, one can expect these results, since Hi's take the same form 
as in general relativity. At this moment, no new constraints are generated, and we have 
to compute the last Poisson bracket {J d 3 x£H, f d 3 yrfH}-ph- This Poisson bracket is too 
complicated to analyze, instead, we consider a simple case that is y/g£, = <5 3 (x — y) with 77 
an arbitrary scalar function. In the process of computation, the following two formula are 
used, 

SRij = ^(ViVjSgki + ViVi§g kj - ViV k 5 9ij - ViVjSgu), (2.22) 
5T) k = ^(VjSga, + V k 5 9jl - V,% fc ). (2.23) 

After a tedious calculation, we obtain the following result expanded with respect to the 
covariant derivatives of rj in different orders. 

{H(x), / d 3 y V H} Ph = -2^-^(a^' fc V fc V i V i r ? + /^V,V,r, + fV^ + urj), (2.24) 
with a l i k given by 

a ijk = ^ilmgjk + cklm g ij _ gilkgjm _ C kli ff> m )K lm , (2.25) 

where C* jfc is defined as e^ l C, k , in which C k = g i m C mk with C mk the Cotton tensor 
defined before. The term including the second order covariant derivative of rj takes the 
form 

FVjVi-n = VyV^Vfcj^^^V* " K lm C ak g n, (2-26) 
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where the notation (kji) c = (ijk + jki + kij) represents the cyclic permutation among 
indices i,j,k. The term proportional to the first order derivative of r/ is a little lengthy, 
and we have to introduce some abbreviations 

j.ijklm _ Qijnigkl _|_ Qivnj gkl _|_ (jjmlgik^ (2.27) 
s ijklm _ jjijrrigkl _|_ fjimj gkl _ fjjml gik ^ 28) 



Then 



f = t mZ ^ i V (i ryV m V fc)c K j7 + ^ ; V(flV,V m)c t ml ^' ! 

+2 S '™^V i r ? V^V fc] ^ m + 2^ m V J 7 ? V [fc V i]S /m ^ 
+2(C klj R\ + <5 /fci i?/ + C lij Rf)K jk Vir] 

H\R)kiC Mj ~ Riu&^KVw (2.29) 



The last one u> is 



u = V l (C^ kl RlK jl + C^ ik R)K kl + C^RlKjt) 

k fx-! _V7 _T7 . i V7T7.V7 f/i V7 V7^ 



1 .i 



-VfcVjVi + K l k VjViVi - K jk V l VtVi - KV 3 V k Vi)C^ k . (2.30) 



At first sight, one may think that constraints from the third covariant derivatives of 77 
are a 2jfc =0, it is not true. The reason is the following. We first notice that 

VkVjVii] = VkViVjT]. (2.31) 

This originates from the formula of second order covariant derivatives of a scalar function, 
which is 

VjViV = djdw - rfjdkri, (2.32) 

where is the Christopher symbol. For a torsionless space, Tfj is symmetric in the 
two lower indices, combining with commutativity of partial derivatives, we have Eq.( 2.3l| ). 
Furthermore, we have the following identities, 

VkVjViV = ^V {k VjV i)c v + g(V fc ViVj - ViV k Vj)v + ^{VkVjVi - VjVfcVi)i7, (2.33) 

where we used Eg. (Hp). With the help of Eq. flOlD , we find that the first term in above 
equation is completely symmetric in three indices. If we denote the symmetrization among 
three indices by (ijk) = ^(ijk+jki+kij+jik+kji+ikj), then |V(fcVjVj) c r? = V(feVjVj)r/. 
The remaining two terms in Eq. (|2.33| ) can be reduced to the first order derivative of r\ by 
utilizing the following formula, 

(V fc ViVj - ViV k Vj)v = Rjik^lV- (2-34) 

After this process, it is clear that the third order covariant derivative term in Eq. (]2.24| ) 
can be written as 

a^VflfeVj-Vi)^. (2.35) 
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Then inheriting the symmetry of V(fcVjVj)ry, the effective components of a % i k are a™ fe '. 
Explicitly, they take the following form 

a (ijk) = 1^klm g ij + (jjim gl k + c ilm gi k )K ml - -(C ilk g jm + C kli g jm + C jli g km )K ml 
3 3 

_1 ^ilj g km + + ci lk g im )K ml . (2.36) 

Because covariant derivatives of r\ of different order are independent, consistency requires 
that the coefficients in front of covariant derivatives of r/ should vanish (Because 7i is 
a first class constraint and the reason will be given later). We deduce from the term 
a^' fc V (fc V j V i) r ? that 

= 0. (2.37) 

To see whether this condition gives rise to new constraints, we will work in a special frame 
where Eq.(^37|) b ecomes simple. First we write the 3-dimensional space metric tensor 
gijdx l dx :i in terms of vielbein, namely 

g ij dx i dx j = 5 ah e a e\ (2.38) 

where 5 a b is the Kronecker delta function, new basis 9 a is related to dx 1 by 6 a = Ofdx 1 . 
Since Cotton tensor is a symmetric tensor, we can diagonalize it through an orthogonal 
transformation O(x) at each given point. So in terms of new basis defined by 6 a = O a b 6 b , 
the component of Cotton tensor becomes 

C afi = {OCO T ) a P = diag{Ci, C 2 , -Ci-C 2 }, (2.39) 

where Ci,C2,and — C\ — C2 are three eigenvalues of Cotton tensor, where we used the 
property that Cotton tensor is traceless. Contracting two sides of Eq. (|2.6| ), we obtain 

VgK = ~ fr = si ii n j ), (2-40) 

this helps us to reexpress Kij by it 13 

VgKV = ^ - X -g^-K. (2.41) 
Now, we find that Eq. flOgp gives the following seven independent equations 

d7r 12 = 0, C 2 7r 12 = 0, (2.42) 

C l7 r 13 = , C 2 tt 13 = 0, (2.43) 

C l7 r 23 = , C 2 n 23 = 0, (2.44) 

(Ci + 2C 2 )ir u - (2d + C 2 )tt 22 + {C x - C 2 )ir 33 = 0. (2.45) 
The solutions of these constraints can be separated into two classes 
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• 1) C\ = C2 = or in other words Cotton tensor should vanish. We notice that the 
other coefficient tensors also vanish because they all consist of Cotton tensor and its 
covariant derivatives. Since the fact that a tensor vanishes is a coordinate independent 
statement, these constraints are still required to be satisfied, even after one uses the 
three coordinate transformations to fix three components of gij. In a special frame, 
we see that the vanishing of Cotton tensor gives two constraint equations, actually, 
in a general coordinate frame, the vanishing of Cotton tensor indeed provides only 
two constraints. The reason is that the following properties of Cotton tensor 

Ca = C iiy g i3 C i3 =0, VjC i3 =0, (2.46) 

make Cotton tensor itself have only two independent components. Gauge-fixing and 
vanishing of Cotton tensor altogether make only one degree of freedom in g^ be 
physical. Now, the Hamiltonian constraint and three momentum constraints elimi- 
nate four conjugate momenta, leaving two components of tt 13 be physical. Altogether, 
the phase space is described by three unpaired fields. 

• 2) 

Ci ^ or C 2 + 0, 7T 12 = 0, vr 13 = 0, vr 23 = 0. (2.47) 

This case is rather bad, it indicates that all the conjugate momenta of g^ are unphys- 
ical, since three momentum constraints already elminates three of the six conjugate 
momenta. 

To interpret above results with the approach given by Dirac , we consider the time derivative 
of the constraint. Utilizing Eqs. (|2T20|) , ( fl2~l"l ) and (plMj) , we obtain 



^ = d l MH + (diNnHj + dj (N j Hi), (2.48) 

^ = (diN^H + AAA, (2.49) 
at 



where the operator A is defined as 

1 



A = -2^57^ (a i3k V fc Vj Vi + jf 3 Vj V< +7^ + 0/), (2.50) 



where the coefficients a l3k , ft 13 , 7* and ui take the same expression as in Eq.(2.24). The 
preservation of constraint in time require ~ 0, ^ ~ 011 the constrained phase 
space. The first one is satisfied because H ~ and TLi ~ 0. While to satisfy the second 
one gives a differential equation AAA ~ 0. It is remarkable that A has no inverse on 
the whole constrained phase space, due to the existence of configuration gij, it 13 making 
the coefficients in front of covariant derivative vanish. As we have analyzed before, the 
configuration g^ with Cij = can achieve this. Therefore, TL cannot be perceived as a 
second class constraint, since to define the Dirac bracket associated with the second class 
constraint, the inverse of A is indispensable. So 7i can only be a first class constraint, and 
its Poisson bracket generates constraints according to the property of first class constraint. 
This justifies our previous treatment of the terms yielded by {Tt(x), f d 3 yr]7i}p\ ) . 
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To summarize, we have found new constraints generated from the Poisson brackets 
of Hamiltonian constraint. These new constraints reduce further the phase space in a 
way that it appears that no symplectic structure exists, or eliminate all the degrees of 
freedom. We expect that the Poisson brackets among the new constraints and 7i, 7ii yield 
more constraints, until the constraints form a closed algebra. When this is done, all the 
constraints are called the first class. Note that the Hamiltonian density is not a second 
class constraint, since the Poisson bracket obviously does not have an inverse. 

Put together, all new constraints either eliminate all degrees of freedom, or make the 
reduced phase space unphysical. 

A complete discussion deserves to be carried out in another work. 

3. Discussion and conclusion 

In previous section, our calculation shows explicitly that Ti and Ti% do not form a closed 
algebra as what happens in general relativity. Intuitively, we feel that non-closure of 
the Poisson brackets among Ti and Ti% can be interpreted by considering the relationship 
between constraints and gauged symmetry . The details are what follows. One can calculate 
the Poisson bracket between the combined constraint NTi + N % Tii and gij then obtains 

{ J d 3 xrj(NH + J^Hi), ^}p b = -divNj - d^N, - V gi,. (3.1) 

It is nothing but the variation of under (3+l)-dimensional coordinate transformation 
in time direction. So the four generators of (3+l)-dimensional coordinate transformation 
are (NTi + N l Hi) and 7ii. This seems to suggest that the constraints Ti and Hi require 
the theory to have full (3+l)-dimensional covariance to respect them, or more constraints 
should be added. The reason is that in a field theory, a local constraint is always accom- 
panied by a gauge symmetry, if we denote the constraint by C, the meaning of constraint 
can be expressed by 

C|phys) « 0. (3.2) 

This equation implies that the physical state is gauge invariant if there is a corresponding 
symmetry. In H-L theory, although there is no diffeomorphism invariance in the time 
direction, Eq.(|3,l|) still looks like such a transformation and imposing this constraint is 
somewhat in conflicts of the starting point 2 . 

What we have discussed previously is based on the UV action of H-L theory. To 
complete our discussion, the full action containing the description of H-L theory both in 
UV and IR region should be taken into account. In the IR region, some operators with 
lower mass dimension will become relevant. The detailed balance principle forces the action 
to take the following form ||] 

S = f dtdPx^gNiiKijK*' - XK 2 ) - ^-C^ 

2 Imposing local Hamiltonian constraint, but lacking the diffeomorphism invariance in the time direction 
may cause strong coupling problem 11 in the IR region of H-L Theory. Before the revised version of our 
paper appears, paper fl9[ also found this independently. 
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k w 4 
+ 4(1 / f 3A) ( 1 ^^ 2 + ^wR - 3A^)}, (3.3) 

For this action to be a deformation of Einstein-Hilbert action in the IR region, it is natural 
to express this action in relativistic coordinates by rescaling t, 

xq = ct, (3-4) 

with the emergent speed of light and effective cosmological constant given by 



At the same time, parameter A should be equal to 1 for the kinetic term taking the same 
form as its counterpart in general relativity. Thus to have a real speed of light, Aw should 
be negative. Setting A = 1, the action becomes 



/ dtd 3 x^gN{(K i:j K^ - K 2 ) - ^d 3 C ij 
J k w 

+^ k R u V 3 R l k - ^R %3 W 



+^Ri-A w R + 3A 2 w )}, (3.6) 

The momentum constraints corresponding to this action remain the same form as in 
Eq,( [2.12|) , the Hamiltonian density contains more terms than in Eq.( [2.1l|) . We denote 
the new Hamiltonian constraint by TC to distinguish it from the previous Hamiltonian 
constraint 

n/V9 = -v ij G ljk iK kl -R + 2A + -^CijC 1 ! - J^e^ k R a V 3 R[ + ^R t3 R 11 
9 % k w 4 



2 

R z , (3.7) 



3 M T>2 



32 

where the speed of light has been set to 1 and then the relation Ay/ = — 8//U 2 has been 
used. It is noticed that the first three terms make up of Hamiltonian constraint in general 
relativity. Calculation of Poisson bracket of Ji is similar to the previous one, and we obtain 
the following result, 

{H(x), [ d 3 y V H} Ph = -2V5r^(« iifc V fc V i V i r ? + ^V 3 V iV + j l V iV + Ct/), 

J «v 

+d i (rig i3 Hi) - gVHidjrj. (3.8) 

The last two terms come from the well known result of Poisson bracket between the Hamil- 
tonian constraint of general relativity. In the previous section, our discussion mainly con- 
centrates on coefficient tensor of the third order covariant derivatives of r], we want to see 
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whether this new coefficient tensor will change our results. Then we compute a l ^ k and 
find that it amounts to replacing the Cotton tensor C 13 in a %3 k by C ij - ^Ri. The 

new contribution —^^-R^ comes from term — ^-e l ^ k Ru'S/jR l k , since it also contains the 

W ■ ■ uk 2 

third order derivative of metric gy. If we denote the tensor C*- 7 y~R %3 by E iJ , then the 

vanishing of 5W fc ) becomes the following seven independent constraints, 

(Si - S 2 )vr 12 = , (E 2 - E 3 )vr 12 = 0, (3.9) 

(Ei - S 2 )tt 13 = , (E 2 - E 3 )vr 13 = 0, (3.10) 

(E 1 -S 2 ) 7 r 23 = 0, (E 2 - E 3 )tt 23 = 0, (3.11) 

(S 2 - S 3 )7r n + (E 3 - E^tt 22 + (E x - E 2 )vr 33 = 0. (3.12) 

To obtain above equations, we also work in basis where Ej,- is diagonalized by Ejj = 
diag{Ei, E 2 , E 3 }. the three eigenvalues are independent because usually E 2 - 7 does not 
satisfy the traceless condition. Despite this difference, our analysis can still be applied to 
this case and our conclusion is unchanged, the same difficulty remains. 

One may choose not to impose the Hamiltonian constraint, this contradicts the require- 
ment that the lapse function is a full-fledged function. Moreover, if one does not impose 
this constraint at the beginning, how can one obtain the usual Hamiltonian constraint in 
Einstein theory in the infrared regime? (This constraint is crucial in going back to the 
Newtonian limit) Finally, we want to remark that even if one can come up with some cure 
of this problem, it will be very difficult to come up with modified theory containing a spin 
2 graviton. 

Note added: When this paper was reviewed, we were informed that another paper 
|1C] holds a different opinion on the constraint structure in Hofava-Lifshitz theory. We 



disagree with the result in [|T^] for the following reason. Most of discussion in that work 
is based on a perturbative method, while this method is not suitable for discussing the 
fundamental degrees of freedom of a theory. For instance, the new constraints obtained 
in our paper are at least of the second order perturbations around Minkowski background 
(adopted in |Io|| ), so at the linear level they do not show up. However, these constraints 
already determine how many degrees of freedom a theory can have before one carries out a 
perturbative calculation. In summary, our point of view is that this perturbative method 
may be useful in solving equations of motion but is invalid for counting the number of 
degrees of freedom of a theory. 
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